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Abstract 

We study a generalization of the Heston model, which consists of two coupled stochastic differ- 
ential equations, one for the stock price and the other one for the volatility. We consider a cubic 
nonlinearity in the first equation and a correlation between the two Wiener processes, which model 
the two white noise sources. This model can be useful to describe the market dynamics char- 
acterized by different regimes corresponding to normal and extreme days. We analyze the effect 
of the noise on the statistical properties of the escape time with reference to the noise enhanced 
stability (NES) phenomenon, that is the noise induced enhancement of the lifetime of a metastable 
state. We observe NES effect in our model with stochastic volatility. We investigate the role of the 
correlation between the two noise sources on the NES effect. 
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I. INTRODUCTION 



The presence of noise in physical systems is a well known phenomenon. It is common 
opinion that noise affecting the djTiamics of a system introduces some degree of instability to 
the system itself, but there is evidence that in some cases the noise can increase the stability 
of a system ^^SSQ]- Noise modeling through the use of the stochastic processes formalism 
has applications that involve many systems including physics, biology, ecology and 
even financial markets 0,0, Q|. The most basic model for financial market is the geometric 
Brownian motion ^ . This model has different drawbacks, it cannot reproduce in fact three 
important stylized facts observed in financial time series: (i) the non Gaussian distribution 

nn 

of returns, (ii) the fat tails j9|, llO[ , and (iii) the stochastic character of volatility, which is 
characterized by long range memory and clustering 

0,0, HQ- 

More complex models 

have been developed to reproduce the dynamics of the volatility. It is worthwhile citing 
the ARCH and GARCH models, where the actual volatility depends on the past 
values of squared return (ARCH) and also on the past values of the volatility (GARCH). 
Another class of models use a system of stochastic equations writing the price as a geometric 
Brownian motion coupled with a non-constant volatility described by a second stochastic 
differential equation. The Heston model uses for the volatility a multiplicative stochastic 
process characterized by mean reversion jiol . IfT^ . Both the models presented so far have 
exponential autocorrelation function so they are not able to reproduce quantitatively the 
long range memory observed in real markets. Nonetheless using values of the characteristic 
time that are sufficiently high, they are able to give accurate statistic for the stock prices, 
by tuning only few parameters. Another important characteristic of financial markets is the 
presence of different regimes. Markets indeed present days of normal activity and extreme 
days where very high or very low price variations can be observed. These are known as 
crash and rally days. A nonlinear Langevin market model has already been proposed 
where different regimes are modelled by means of an effective potential for price returns. In 
some circumstances this potential has a cubic shape with a metastable state and a potential 
barrier. The dynamics inside the metastable state represents the days of normal evolution 
while the escape after the potential barrier represents the beginning of a crisis. Metastable 
states are ubiquitous in physics and the effect of noise in such systems has been extensively 
studied , but considering the noise intensity as a parameter (physical models 
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described by additive stochastic differential equations for example). Financial markets with 
their stochastic volatility are an example of systems where the noise intensity is far from 
being a constant parameter, but it is indeed a stochastic process itself. Moreover there is 



evidence in nature that t 
multiplicative noise 



le noise intensity is not a constant parameter and can be modelled 



19|. So it is interesting to release the hypothesis of parametric 



noise intensity in financial market models as well as in natural systems. 

II. THE HESTON MODEL WITH A METASTABLE STATE 

The Heston model introduced in the previous section is described by the following system 
of coupled stochastic differential equations [3| 



dx{t) = {ii-v{t)/2) ■ dt + ^/v{f) ■ dZ{t) 



dv{t) = a{b-v{t))-dt + c^/v{t)-dW{t) 



(1) 



The price p{t) follows a geometric random walk whose standard deviation is another 
stochastic process. Here x{t) = Inp(t) is the log of the price, Z{t) and W{t) are un- 
correlated Wiener processes with the usual statistical properties: (i) {dZ{t)) = and 
{dZ{t) ■ dZ{t')) = 6{t - t')dt; (ii) {dW{t)) = and {dW{t) ■ dW{t')) = 6{t - t')dt. The 
v{t) process is characterized by mean reversion, i.e. its deterministic solution has an ex- 
ponential transient with characteristic time equal to a^^, after which the process tends to 
its asymptotic value b. The process for v{t) exhibits the phenomenon of volatility cluster- 
ing, alternating calm with burst periods of volatility, and has an exponential autocorrelation 



function. The smaller the value of a the longer are the bursts in vo 
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atility. Heston model has 
23[ . The equations of the 



been subject of recent investigation by econophysicists 
system (1) are well known in finance, they represent respectively the log-norm al g eometric 
Brownian motion stock process used by Black and Scholes for option pricing j^, H]; ^^'^ 
the Cox-Ingersoll-Ross (CIR ) rn ean-reverting stochastic differential equation first introduced 



for interest rate models 
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Here we consider a generalization of the Heston model, by replacing the geometric Brow- 
nian motion with a random walk in the presence of a cubic nonlinearity. This generalization 
represents a "Brownian particle" moving in an effective potential with a metastable state. 
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in order to model those systems with two different dynamical regimes like financial markets 
in normal activity and extreme days Q]. The equations of the new model are 



dx{t) 



dv{t) = a{b-v{t))■dt + c^/vit)■dW{t), (2) 
where U is the effective cubic potential U (x) = px^ + qx"^, with p = 2 and q = 3 (see Fig. HI), 



Z{t) and W{t) arc 



"if the potential 




FIG. 1: Cubic potential used in the dynamical equation for the price x(t). The points in the figure 
indicate the starting positions used in our simulations. 

maximum and xi the abscissa where the potential intersects the x axes. The intervals 
Xo < xj and / = [xj, xm] are clearly regions of instability for the system. In systems with a 
metastable state like this, the noise can originate interesting effects that increase instead of 
decrease the stability, by enhancing the lifetime of the metastable state [l|,flS3. Agood 
example is the Noise Enhanced Stability (NFS) phenomenon. The mean escape time r for 
a Brownian particle moving throughout a barrier AU, with a noise intensity v, is given by 
the well known exponential Kramers law 13, Isol 



Aexp 



AU 



(3) 



where r is a monotonically decreasing function of the noise intensity v, and A is a prefactor 
which depends on the potential profile. This is true only if the random walk starts from 
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30|. For a detailed discussion on this point and different dynamical regimes 



initial positions inside the potential well. When the starting position is chosen in the insta- 
bility region Xo < xm, t exhibits an enhancement behavior, with respect to the deterministic 
escape time, as a function of v. Particularly for initial positions Xq < xj, we have nonmono- 
tonic behavior of r as a function of v 0, . This is the NES effect and can be explained 
considering the barrier "seen" by the Brownian particle starting at the initial position Xq, 
that is AUin = U{xmax) — U{xo). Moreover AUin is less than AU as long as the starting 
position Xq lyes into the interval I = [xi,xm]- Therefore for a Brownian particle, from a 
probabilistic point of view, it is easier to enter into the well than to escape from, when the 
particle is entered. So a small amount of noise can increase the lifetime of the metastable 

state iy,.^^ 
TlTT 

see Refs. mm. When V is much greater than AU, the Kramers behavior is recovered. The 
NES effect has been experimentally observed in a tunnel diode and theoretically predicted 
in a wide variety of systems such as for example chaotic map, Josephson junctions, chemical 
reaction kinetics, and neuronal dynamics models . Our modified Heston model, 

characterized by a stochastic volatility and a nonlinear Langevin equation for the returns, 
has two limit regimes, corresponding to the cases a = 0, with only the noise term in the 
equation for the volatility v{t), and c = with only the reverting term in the same equation. 
This last case corresponds to the usual parametric constant volatility case. In fact, apart 
from an exponential transient, the volatility reaches the asymptotic value b. The NES effect 
should be observable in the latter CclSG clS 8b function of b, which is the average volatility. In 
this case, in fact, we have the motion of a Brownian particle in a fixed cubic potential with 
a metastable state and an enhancement of its lifetime for particular initial conditions (see 
Refs. BQ). 

III. ENHANCEMENT OF THE ESCAPE TIME 

The two processes of Eqs. (JJ) and are actually uncorrelated. In financial markets the 



two processes can 
as leverage effect 
scientific literature 



3e correlated, and a negative correlation between the processes is known 
31l. Heston model with correlation has been recently discussed in the 
Our modified Heston model becomes therefore 
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dx{t) = -(^ + -Lij .dt + ^/^)-dZ{t) 




dv{t) = a{h - v{t)) ■ dt + c^/v{¥) ■ dWcit) 



dWc{t) = p-dZ{t) + -dWit), 



(4) 



where Z(t) and W(t) are uncorrelated Wiener processes as in Eqs. ((T)) and (j21), and p is 
the cross correlation coefficient between the noise sources. The investigation is performed 
simulating the process of Eqs. Q with time integration step At = 0.01, and for the fixed 
starting position xq = —1.25 in the I interval (this initial position is shown as a black square 
point in Fig. The absorbing barrier is located at x = —6.0, and the results are averaged 
over 10^ escape events. The algorithm used to simulate the noise sources in Eqs. (0} is based 
on the Box-Muller method for generating random processes with a Gaussian distribution. 
The numerical integration of Eqs. Q is done by using the forward Euler method js^. 

Our first result shows that the curve r vs. h is weakly dependent on the value of the p 
parameter. This is shown in Fig. |21 where all the curves correspond to the region of the 
parameters space where the effect is observable. There is indeed a weak variation in the 
maximum value. The highest maximum values correspond to the highest absolute p values. 
It is worth noting that the correlation affects directly the noise term in the v{t) equation, 
and it has negligible influence on the reverting term of the same equation. 

The curves of r vs. c however have an evident dependence on the parameter c. This is 
shown in Fig. El for two values of parameter a for which the NES effect is observable. The 
correlation affects the position as well as the value of the maximum of r (see Fig. IH}, but in 
a different way. The position of the maximum increases only for very high positive values of 
the correlation coefficient. The maximum value of r increases with positive correlation but 
decreases with negative correlation. This effect is more evident for the higher value of a, as 
shown in Fig. Et". To comment this result we ffist note that a negative correlation between 
the logarithm of the price and the volatility means that a decrease in x{t) induces an increase 
in the volatility f (t), and this causes the Brownian particle to escape easily from the well. As 
a consequence the mean lifetime of the metastable state decreases, even if the nonmonotonic 
behavior is still observable. On the contrary, when the correlation p is positive, the Brownian 
particle stays more inside the well, decrease in x{t) indeed is associated with decrease in the 
volatility. The escape process becomes slow and this increases further the lifetime of the 
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FIG. 2: Mean escape time r as a function of the mean volatility b. The other parameters have 
the fixed values a = 10~^ and c = 10~^. The six curves correspond to the following values of 
the correlation: p = —0.8, —0.5, —0.1,0.1,0.5,0.8. The fixed starting position is xq = —1.25. The 
values of the potential parameters (see Eq. (j2j) are: p = 2, q = 3. 




Q G 

FIG. 3: Mean escape time r as a function of the model noise intensity c for two different values 
of the parameter a: (a) a = 2.0 and (b) a = 20. The value of the parameter b is fixed to b = 10~^. 
The fixed starting position xq and the potential parameters p and q are the same of Fig. |2| The 
different curves correspond to the following values of p: -0.8 (circle), 0.0 (square), 0.8 (triangle). 
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metastable state, causing an increase in the value of the maximum of the curve of Fig. 01 
To illustrate better this aspect and the behaviour of Figs. |21 and El near the maximum, we 
plot the values Tmax of the maximum as a function of the correlation coefficient p and we 
show these curves in Fig. |31 Specifically in Fig. we report the values of T^ax related to 
the curves of r of Fig. |21 and in Fig. ^ those related to Fig. |21 The increase of stability is 



' max 




FIG. 4: The values of 

Tmax & function of the correlation coefficient of the curves reported in 
Figs. 121 and O Specifically: (a) Tmax from Fig. (21 (r as a function of parameter b) with a = 10~^, 
c = 10^^); (b) Tmax from Fig. |21 with a = 2.0 (circles) and a = 20 (squares), b = 10~^. The fixed 
starting position xq and the potential parameters p and q are the same of Fig. [21 



evident in panel (b) as Tmax increases with p. 

It is interesting to show the probability density function (PDF) of stock price returns for 
the model described by the Eqs. (El). This is done in Fig. El As one can see the qualitative 

nn 

behavior of a fat tail distribution, typical of real financial market data, is recovered [21, [10| , 
but with a peculiar asymmetry. To characterize quantitatively the PDF of returns (Fig. \^ 
as regards the width, the asymmetry and the fatness of the distribution, we calculate the 
mean value < Aa; >, the variance a^x, the skewness and the kurtosis K4. We obtained the 
following values: < Ax >= -0.434, aAx = 0.903, K3 = -2.086, K4 = 9.417. These statistical 
quantities clearly show the asymmetry of the distribution and its leptokurtic nature observed 



in empirical investigations, characterized 



tails than in the Gaussian distribution 



jy a narrow and larger maximum, and by fatter 



The presence of the asymmetry is very interesting and it will be subject of future inves- 
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FIG. 5: Probability density function of stock price returns for the Heston model with a metastable 
state and correlated noise sources (Eqs. @). The parameters of the Cox- IngersoU- Ross process of 
Eqs. (0J) are: a = 2, b = 0.01, c = 26. The cross correlation coefficient is p = —0.3. The potential 
parameters are the same of Fig. |21 



tigations. It is worthwhile to note, however, that the PDF of returns become asymmetric 
in crash and rally days that are just the time periods of financial data where the 

related dynamical regimes could be described by the models (j2I) and (jH). Of course the 
quantitative agreement between the PDF of real data and that obtained from these models 
requires further investigations on the parameter value choice of the cubic potential, the pa- 
rameters of the simple Heston model and the correlation coefficient. This analysis is outside 
the aim of the present work and will be further investigated in a forthcoming paper. 

Finally in the following Fig. El we report the comparison between the probability density 
function of the escape times of daily price returns from real market data and that obtained 
from the model described by the Eqs. (jH). The data set used here consists of daily closure 
prices for 1071 stocks traded at the NYSE and continuously present in the 12-year period 
1987 — 1998 (3030 tradin g d ays). The same data were used in previous investigations by 
one of the authors [2^ 0,]^ . From this data set we obtained the time series of the returns 
and we calculated the time to hit a fixed threshold starting from a fixed initial position. 
The two thresholds were chosen as a fraction of the standard deviation cr^ observed for 
each stock during the above mentioned whole time period (n is the stock index, varying 
between 1 and 1071). Specifically we chose: (Axj)„ = — 0.1cr„ and (Ax/)„ = — 1.0(T„. The 
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parameters of the CIR process are: a = 20, b = 0.01, c = 2.4. The cross correlation 
coefficient is p = —0.9, and the potential parameters are the same of Fig. |21 As one can 
see the agreement between real and theoretical data is very good, except at small escape 
times. The choice of this parameter data set is not based on a fitting procedure as that used 
for example in Ref. ||21|, where the minimization of the mean square deviation between the 
PDF of the returns extracted from financial data and that obtained theoretically is done. 
We chose the parameter set in the range in which we observe the nonmonotonic behaviour 
of the mean escape time as a function of the parameters b and c. Then by a trial and error 
procedure we selected the values of the parameters a, b, and c for which we obtain the best 
fitting between the PDF of escape times of the price returns calculated from the modified 
Heston model (Eqs. (j^) and that obtained from time series of real market data. Of course a 
better quantitative procedure could be done, by considering also the potential parameters. 
This will be done, together with a detailed analysis of PDF of returns and its asymmetry, 
in a forthcoming pa^^'^'' 

10-' 
10-' 

Q 

'^lo-^ 
10-^ 

50 100 

escape time (days) 

FIG. 6: Comparison between the probability density function of escape times of the price returns 
for the modified Heston model with a metastable state (Eqs. (jlj)) (soHd line) and the PDF of escape 
times of returns for real market data (cirles). The parameter values of the model @ are: p = 2, 
q = 3, a = 20, b = 10~^, and c = 2.4. The cross correlation coefficient between the noise sources 
is: p = -0.9. 



» real data 
— theoretical results 
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IV. CONCLUSIONS 



We have investigated the statistical properties of the escape time in a generahzed Heston 
market model, characterized by the presence of a metastable state in the effective potential 
of the logarithm of the price p{t). We observe the NES effect in the system investigated. 
The presence of correlation between the stochastic volatility and the noise source which 
affects directly the dynamics of the quantity x{t) = Inp(t) (as in usual market models), can 
influence the stability of the market. Specifically a positive correlation between x{t) and 
volatility v{t) slows down the walker escape process, that is it delays the crash phenomenon 
by increasing the stability of the market. A negative correlation on the contrary accelerates 
the escape process, lowering the stabihty of the system. 

This work was partially supported by MIUR. 
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